Plant Direct

Water-use efficiency (WUE), the ratio of accumulated plant biomass to water lost through transpiration has conventionally been determined using a destructive single-point measurement. Recent advances in high-throughput phenotyping now enable repeated, non-destructive estimation of biomass and WUE. However, these digital measurements must be statistically validated against conventional destructive methods to validate their use as reliable proxies. Therefore, we compared digital biomass determined point clouds produced from multispectral camera scanners with destructive harvests across eight harvests using Samsun tobacco grown under both drought and high-water conditions. WUE efficiency, calculated using the digital biomass estimated from a point cloud and gravimetric water use determinations, were compared to destructive harvest determinations. The coefficient of variation (CV) showed there were no significant differences in digital and destructive measurements for either biomass or WUE. Indicating that digital measurements can be used in place of destructive measurements. Drought plants used significantly less water and were significantly smaller than high-water plants from Harvests 4 through 8. However, there were no significant differences in the ratio of evapotranspiration to leaf area or WUE, indicating that drought plants were simply smaller and used less water than the high-water plants. This work validates that estimating plant biomass from a digital point coupled with continuous gravimetric determination of water use provides a reliable nondestructive measure of WUE in high-throughput measurements across the full plant life cycle.

High-phenolic plant genera present a major technical limitation in genomic research. Standard extraction approaches that perform reliably across diverse flora often perform poorly when applied to recalcitrant taxa, producing low DNA yield and integrity incompatible with sequencing requirements. The genus Anigozanthos (Kangaroo paws) from the family Haemodoraceae exemplifies this problem. We identified key physicochemical factors governing extraction failure in this genus and resolved them through targeted modifications to lysis chemistry and contaminant management. The resulting protocol achieved a near threefold improvement in DNA purity, substantially reducing contaminant carry over and consistently yielded high-integrity, long DNA fragments (DIN > 7) across a diverse sample set spanning cultivated and wild material across four diverse genera of Haemodoraceae. We also tested a straightforward purity assessment framework that can be implemented in any standard molecular laboratory, enabling rapid pre-submission quality assessment without the need for specialised equipment. Together these advances open a practical path to genomic characterisation of Anigozanthos that establishes a transferable model for genomic research across Australia ' s chemically complex native flora.

Although the biosynthetic pathways of benzylisoquinoline alkaloids (BIAs) have been extensively investigated in several plant species, their transcriptional regulatory mechanisms remain only partially understood. Jasmonate (JA)-responsive group IX APETALA2/Ethylene Responsive Factor (AP2/ERF) transcription factors (TFs) are well-known regulators of specialized plant metabolism, including the biosynthesis of various alkaloids. However, their specific roles in BIA biosynthesis remain largely elusive. Here, we isolated five novel group IX AP2/ERF TFs, designated Benzylisoquinoline alkaloid Jasmonate-responsive AP2/ERF (BJE1-5), from Coptis japonica. Phylogenetic analysis revealed that Benzylisoquinoline alkaloid Jasmonate-responsive AP2/ERF (BJE) proteins belong to subclades distinct from group IXa, which contains well-known AP2/ERF TFs involved in alkaloid biosynthesis. Transient expression analyses in C. japonica protoplasts demonstrated that certain BJEs, particularly CjBJE3 and CjBJE5, positively regulated BIA biosynthetic genes through a mutual regulatory network among BJE members. Moreover, CjBJE3 expression was regulated by CjbHLH1, a unique-type basic helix-loop-helix (bHLH) TF specific to BIA-producing plants. Furthermore, heterologous expression of CjBJE3 and CjBJE5 in cultured Eschscholzia californica cells significantly enhanced the overall BIA production, particularly by increasing end-product benzophenanthridine BIAs, highlighting several uncharacterized biosynthetic genes clustered in the genome. Our findings suggest that BIA-producing species have developed a specific regulatory network comprised of CjbHLH1 and BJE TFs, providing valuable clues for identifying novel biosynthetic enzymes.

MicroRNAs (miRNAs) are key post-transcriptional regulators of development, stress response, and secondary cell wall formation in woody plants, yet annotations for Eucalyptus grandis, the world's most widely planted hardwood, remain fragmented across studies using incompatible discovery pipelines and filtering criteria. Here we present the Eucalyptus MicroRNA Archive (EMA), a curated, locus-resolved database integrating three independent small RNA sequencing datasets spanning vegetative tissue, somatic embryogenesis, and mechanically induced tension wood formation. Applying annotation criteria aligned with current plant miRNA standards, EMA catalogs 99 curated miRNAs (31 previously described, 68 novel) organized into 34 family-level groupings under a three-tier confidence system, known-reference-supported, multi-study replicated, or single-study, that preserves study-of-origin and sample-level evidence for every entry. Cross-study comparison showed that only 9 of 99 entries (9.1%) were independently supported by all three datasets, supporting an evidence-tiered rather than binary annotation scheme. Target prediction against the E. grandis transcriptome yielded 1,773 miRNA-target interactions spanning 764 loci, integrated into a combined miRNA-target and protein-protein interaction network. This network resolved into functionally coherent, mutually isolated clusters, including an miR482-associated NBS-LRR/TIR disease-resistance hub with a substantial translational-repression component, alongside modules enriched for ribosome biogenesis and translation, DNA replication, and nitrogen and carbohydrate metabolism. EMA is publicly accessible through an interactive web dashboard, with all curated data, source code, and analysis scripts openly available, providing a reproducible, extensible framework for E. grandis miRNA research and a template for similarly structured resources in other non-model woody species.

Accelerating photoprotective regulation to improve carbon assimilation is a promising strategy to increase crop productivity. Although rapid non-photochemical quenching (NPQ) relaxation has been validated as a target through metabolic engineering, it remains unclear whether conventional breeding has improved this trait. Here, we investigated whether more than a century of soybean breeding enhanced NPQ relaxation alongside light-saturated carbon assimilation and seed traits. We evaluated a historical panel of 24 soybean genotypes across vegetative and reproductive developmental stages by integrating NPQ relaxation, gas exchange parameters, xanthophyll-cycle pigment profiles, expression of key photoprotective genes (VDE, PsbS, and ZEP), seed number and seed weight. NPQ relaxation parameters were not consistently associated with genotype release year, seed number, or seed weight at either developmental stage. The only exception was the amplitude of the rapidly relaxing NPQ component (AqE), which was negatively correlated with all three variables during the reproductive stage. In contrast, genotype release year was positively associated with maximum net CO2 assimilation rate (Amax), maximum carboxylation rate of Rubisco (Vcmax), maximum electron transport rate (Jmax), seed number, and seed weight, while Amax and Vcmax were positively correlated with seed number and seed weight. These findings indicate that the greater photosynthetic capacity of modern genotypes was not accompanied by faster photoprotective response. Thus, photoprotective regulation has not kept pace with gains in photosynthetic capacity under field conditions. We conclude that rapid NPQ relaxation remains an important target for synchronizing photoprotection with the high photosynthetic capacity of modern soybean lines.

All living organisms rely on the movement of ions across cell membranes as the fundamental physical basis of their internal energy and signaling, and plants are no exception. Plants perceive, integrate, and respond to environmental stimuli through electrical signals, classified as action, variation, and system potentials, that are coupled with calcium waves, reactive oxygen species, and hydraulic and hormonal changes to coordinate whole-organism responses despite the absence of a nervous system. Yet most studies characterize these signals using a single feature, such as amplitude or spike duration, in a single tissue, an approach that cannot establish how such signals correspond to the underlying ionic activity, mobility, and structural complexity of the signaling environment, or how this correspondence varies across organs. Here, we correlate plant bioelectrical signals with potential ionic energy flow using a multi-domain framework, combining discrete spike events, continuous waveform properties, spectral composition, and signal complexity applied to leaf, stem, and root recordings from tomato (Solanum lycopersicum) exposed to different stimulus. Electrical activity with increased stimulus strength, likely reflecting increased ionic flow, with the root showing the largest response. This suggests plant electrical signaling works as a distributed, ion-based information system, useful for stress monitoring and bio-inspired sensor design.

Perennial grains represent a promising frontier for sustainable agriculture, but breeding progress is constrained by the accessibility of genotyping and the difficulty of evaluating complex traits expressed for multiple years after establishment across heterogeneous environments. Phenomic selection may help address these challenges by using inexpensive, scalable, high-dimensional phenotypes collected early in development, although the robustness of such predictions across breeding cycles remains uncertain. Here, we compared genomic selection and phenomic selection across two breeding cycles of Thinopyrum intermedium (intermediate wheatgrass; IWG; Kernza(R)), comprising approximately 2,280 individuals from maternal half-sib families evaluated across multiple field sites and years. We constructed relationship matrices from genomic markers and early-life stage phenomic data, including seed and leaf color (HSV), CropReporter multispectral reflectance and indices, and cycle-specific hyperspectral reflectance sensors. Genomic models provided the strongest predictions on average across all field traits in both cycles. Among phenomic predictors, leaf HSV was consistently the most informative, whereas CropReporter and hyperspectral data showed lower and more trait-dependent performance and seed HSV provided little predictive value. Genomic, leaf HSV, and CropReporter models transferred across breeding cycles with little apparent loss of predictive ability relative to within-cycle validation, demonstrating that their predictive signals were not restricted to a single breeding cycle. Early-life stage leaf HSV emerged as a practical, accessible tool for germplasm thinning and early-stage prioritization in perennial breeding programs. Despite limited similarity among relationship matrices, multi-relationship-matrix models rarely improved prediction beyond the stronger constituent single-relationship-matrix model. Together, these results show that early-life stage phenomic data provide reproducible information about agronomic performance expressed years later, but that predictor complexity and data integration do not guarantee improved prediction.

Three new species of sun orchid (Thelymitra) endemic to Aotearoa New Zealand are here described. These are T. palustris, T. scabrifolia and T. semaphora. The morphological distinctiveness of these three species has been acknowledged for decades; however, their taxonomic status has remained unresolved. Evidence from existing karyological data, recently generated DNA sequence data (LFY and ycf1) and morphological studies from historical and fresh collections are used here to support their formal description. Both, T. palustris and T. semaphora are restricted to wet habitats north of Auckland (North Island). Thelymitra scabrifolia inhabits mostly scrub, and it has a similar northern North Island distribution, but is has been found also in Manawat[a]whi / Three Kings Islands and historically in Otago (South Island). All three species are polyploids and are of conservation concern.

Plant cold acclimation emerges from coordinated adjustments in photosynthesis, primary metabolism, and intracellular carbon allocation. Yet, the regulatory role of subcellular metabolite compartmentation in natural variation of cold acclimation remains insufficiently understood. Here, we investigated four Arabidopsis thaliana accessions grown either individually or in bulk to determine how growth configuration and genotype shape the metabolism of sugars and organic acids during cold exposure. Using non-aqueous fractionation, we quantified plastidial, cytosolic, and vacuolar sugar pools alongside whole-cell carbohydrates, organic acids, enzyme activities, photosynthetic parameters, and stress markers. A neural-network classifier revealed that subcellular sugar distribution together with sugar amounts and organic acids provided the strongest discriminatory power among accessions, surpassing photosynthetic traits and enzyme activities. Our findings demonstrate that natural variation in cold acclimation is strongly determined by genotype-specific subcellular metabolite architectures, and that the cultivation strategy modulates these intracellular signatures. We conclude that subcellular compartmentation of metabolites represents a cellular control layer for natural variation of cold acclimation and resilience in Arabidopsis thaliana.

Extrafloral nectaries (EFNs) are a widespread plant defense mutualism trait and are highly convergent, appearing in hundreds of plant lineages worldwide. Here we investigate a report of possible EFNs in Erythranthe angulosa, a recently-described California wildflower. We integrate field observations, insect bioassays, an induction experiment, and microscopy to test for signatures of EFN function, finding no evidence that the distinctive axillary swellings produced by E. angulosa function as EFNs. We also uncovered two distinct morphs at the type locality of E. angulosa that diverge in the number of axillary swellings produced, as well as other shoot architecture traits such as stem thickness, leaf size, and branch number. Although the axillary swellings appear to not function as EFNs, they remain a compelling morphological variant within the yellow monkeyflowers that may perform storage or another unknown function.

BACH1 is a heme-regulated basic-leucine-zipper containing transcriptional repressor that binds its DNA recognition elements as a heterodimer with MAFK. Heme-binding is thought to be mediated by several Cys-Proline (CP) motifs and this results in dissociation of the heterodimer from DNA. The mechanism of heme-binding and heme-mediated DNA dissociation remains unresolved. We have used UV-visible spectroscopy, 2D-NMR and DNA-binding assays to explore both heme-binding and DNA dissociation of a minimal BACH1 construct containing 2 CP motifs (C492(CP5) and C646(CP6)) flanking the DNA-binding domain. We find that heme is able to bind to both CP motifs, but also to other non-CP cysteines and histidines in the construct. Using NMR spectroscopy, we identify a structured binding pocket in which heme interacts with both C646(CP6) and Cys621. However, DNA-binding assays show that C646(CP6) is not required for heme-mediated DNA dissociation of the BACH1:MAFK heterodimer. Using UV-visible spectroscopy we show that C492(CP5) also recruits heme with a second ligand, a conserved histidine, His559, in the BACH1 DNA-recognition helix. Mutation of C492(CP5) reduces but does not abolish heme-mediated dissociation from DNA. Our findings suggest a mechanism for heme-binding to BACH1 and heme-mediated dissociation from DNA.

Anterograde and retrograde signalling establish bidirectional communication between the nucleus and chloroplasts. Retrograde signals from chloroplasts regulate nuclear gene expression while anterograde signals from the nucleus coordinate chloroplast development and maintain cellular homeostasis. How this bidirectional signalling framework extends beyond locus-specific regulation to shape the global spatial organization of nuclear chromatin across tissues remains unclear. Although anaesthesia can alter chromatin organisation, the role of chloroplast dysfunction in these changes remains unclear. Here, we investigate how chloroplast dysfunction and anaesthesia influence euchromatin and heterochromatin organisation in Solanum lycopersicum seedlings across tissues with contrasting photosynthetic competence. Using confocal and super-resolution radial fluctuation (SRRF) imaging with quantitative multiparameter analysis, we identify distinct, tissue-specific chromatin responses to chloroplast disruption and anaesthesia. Notably, anaesthesia induces distinct spatial chromatin changes across tissues that are independent of chloroplast dysfunction, suggesting a direct nuclear response to anaesthesia rather than a chloroplast-mediated retrograde effect. These findings highlight chromatin topology as a potential quantitative biomarker of cellular disruption and provide a framework for investigating anterograde chloroplast-nucleus coordination and stress-responsive nuclear organisation in plants.

Blumeria graminis f. sp. avenae (Bga), the causal agent of oat powdery mildew, is one of the most host-specialized members of the B. graminis species complex. Despite its agricultural importance, the lack of a high-quality reference genome has limited studies of host specialization, virulence evolution and comparative genomics in this pathogen. Here, we generated the first chromosome-scale genome assembly of Bga using an integrative approach combining long- and short-read sequencing, Hi-C scaffolding and transcriptome data. The Bga genome exhibits hallmark features of powdery mildew fungi, including extensive repeat content and low gene density. Comparative analyses revealed that genome expansion is primarily associated with historical transposable element proliferation rather than recent transpositional activity. Genome organization is consistent with a functionally stratified "one-speed" model, in which genes associated with pathogenicity, including predicted effectors and infection-responsive genes, are preferentially located in transposable element-rich regions characterized by reduced synteny conservation and extended intergenic spaces. In contrast, conserved genes are concentrated in compact genomic regions and maintain strong syntenic conservation across cereal-infecting formae speciales. Hi-C analyses demonstrated a highly structured chromatin architecture and revealed genome organization patterns associated with infection-related gene expression. Comparative genomic analyses indicated that host specialization in Bga is driven by localized diversification of a relatively small subset of genes rather than large-scale genome restructuring. These results provide the first high-quality genomic resource for Bga and offer new insights into the evolutionary mechanisms underlying host specialization in powdery mildew fungi.

The integration of digital technologies for high-throughput field phenotyping is critical for accelerating crop improvement in agriculture. However, extracting traits from remote sensing data remains constrained by fragmented workflows, manual intervention, and limited interoperability among existing tools, resulting in delays that hinder timely biological insight and decision-making. To address these challenges, we present PhenoStream (Phenotyping Streaming), a scalable, end-to-end cyberinfrastructure designed to automate the full lifecycle of aerial imagery-based phenotyping, from data acquisition to plot- and genotype-level inference. The framework integrates automated data ingestion from distributed field sites, geospatial processing, and AI-enabled trait extraction within a unified, user-accessible graphical interface. Its modular and extensible architecture supports adaptable trait modeling and seamless integration of new data sources, enabling deployment across diverse crops, environments, and experimental designs. We demonstrate the system across a large multi-location field trial network of bioenergy crops, where it enables high-throughput characterization of spatiotemporal growth dynamics, genotype-by-environment (GxE) interactions, and predictive modeling of key agronomic traits. By significantly reducing processing latency and manual effort, the platform facilitates near-real-time analysis and reproducible workflows. This work establishes a generalizable and scalable pathway for operationalizing very-high-spatial resolution aerial phenotyping in agricultural research. By bridging data acquisition and analytics, the end-to-end cyberinfrastructure provides a foundation for integrating heterogeneous and unstructured data streams--including remote sensing, environmental, and management data--toward data-driven decision making in agriculture.

THAP9 is a transposable element-derived gene which encodes a protein that is homologous to the active Drosophila P-element transposase (DmTNP). Both THAP9 and DmTNP possess a C-terminal domain (CTD) which is functionally uncharacterized. Sequence and structural analysis suggest that the THAP9-CTD has a novel fold which is only found in THAP9 homologs. To explore the evolutionary history and characteristics of this novel domain, exhaustive phylogenetic analysis (using MSA, structure prediction, MSTA-based clustering) was performed. THAP9-CTD homologs were more widely distributed throughout the animal kingdom in comparison to DmTNP-CTD homologs which were restricted to arthropods. Moreover, the THAP9-CTD homologs were more conserved, especially among mammals and birds and their average length increased in a class-specific manner. Comparison with the DmTNP-CTD homologs demonstrates that although their respective CTDs may have evolved independently, they both surprisingly share similar secondary structure elements consisting of three conserved helical regions made of hydrophobic residues that are predicted to make up a conserved core. The role of the respective CTDs were further investigated by creating truncation mutants lacking the CTD. Interestingly both THAP9 and DmTNP truncation mutants are still capable of DNA excision and integration suggesting that their respective CTDs are not essential for DNA transposition. Moreover, CTD truncation favours DNA integration in THAP9: this suggests that CTD acquisition during evolution may have led to THAP9 domestication as observed in other transposable element-derived genes like Rag1 and piggybac, which have similar terminal regulatory domains.

Aegilops peregrina is a wild allotetraploid wheat wild relative and an important source of genetic diversity for stress tolerance and agronomic traits. Here, we report a subgenome-resolved, chromosome-scale reference genome assembly of a drought tolerant and stem rust resistant Ae. peregrina accession PI 604178 generated using PacBio HiFi and Hi-C sequencing. The 10.13 Gb assembly contains 98.81% of sequence anchored to 14 pseudomolecules representing the seven S and seven U chromosomes, with contig and scaffold N50 values of 25.84 and 746.48 Mb, respectively. The assembly achieved a consensus quality value of 74.61, 97.83% k-mers completeness, and 99.9% BUSCO completeness. LTR Assembly Index values of 20.43 and 18.79 for the S and U subgenomes, respectively, further supported high continuity across repeat-rich regions. Repetitive elements comprise 85.93% of chromosome-anchored assembly. We annotated 59,910 high-confidence protein-coding genes, with comparable gene representation across the two subgenomes. This reference genome provides a high-quality genomic framework for comparative analyses, characterization of important loci regulating agronomic and resilience related traits, and sequence-guided exploitation of Ae. peregrina allelic diversity for wheat improvement.

Background: Acute phase response is an early immunometabolic response to brain injuries, primarily coordinated by the liver via the activation of acute phase proteins. These immune responses can be both beneficial, promoting tissue repair, and detrimental, exacerbating neurological deficits, if not properly controlled. Despite the central role of the liver in immunometabolism, how hepatic metabolism dynamically adapts to traumatic brain injury remains under explored, primarily due to limited liver-specific modalities that can assess metabolic pathways in vivo. 13C MRI utilizing hyperpolarized 13C-pyruvate can assess key regulatory enzyme activities in hepatic metabolism. Methods: Rats with controlled cortical impact were studied in vivo using hyperpolarized [1-13C]pyruvate and [2-13C]pyruvate under fed and fasted conditions 3-4 days after injury. Hyperpolarized 13C products, including [13C]bicarbonate from [1-13C]pyruvate and [5-13C]glutamate, [1-13C]acetyl-L-carnitine, and [2-13C]phosphoenolpyruvate from [2-13C]pyruvate, were evaluated to assess mitochondrial and gluconeogenic metabolism. In parallel, liver tissues were collected following [U-13C3]pyruvate injection for NMR isotopomer analysis of phosphoenolpyruvate, glucose, and glutamate. Results: While no metabolic differences were detected under fed condition, [13C]bicarbonate and [2-13C]phosphoenolpyruvate increased after brain injury under fasted condition, indicating an upregulation of the hepatic gluconeogenic pathway after injury. 13C NMR of liver tissue extracts from injured rats showed an elevated [2,3-13C2]glutamate-to-[4,5-13C2]glutamate ratio and increased 13C-labeling in phosphoenolpyruvate than controls, confirming enhanced hepatic gluconeogenic pathway. Conclusion: This study demonstrates that hepatic acute phase response to brain injuries can be monitored in vivo by hyperpolarized pyruvate, which may be further utilized for longitudinal immunometabolic evaluation of the liver during pathogenesis and therapeutic interventions.

The binding affinities of a chimeric analog of a riboflavin derivative linked to biotin, (6- (7,8-dimethyl-2,4-dioxo-3,4-dihydrobenzo[g]pteridin-10(2H)-yl)hexyl 5-((3aS,4S,6aR)-2- oxohexahydro-1H-thieno[3,4-d]imidazol-4-yl)pentanoate), referred to as C6-Rf-biotin-tag, to the riboflavin binding retain or streptavidin are in the M range, 1.29 {+/-} 0.277 and 3.00 {+/-} 0.459, respectively. These values suggest that C6-Rf-biotin-tag has potential applications in diagnostic assay and labelling target flavin binding proteins. The C6-Rf-biotin-tag which was characterized with respect to physical and biochemical properties retains UV/Vis spectroscopic and fluorescence behavior similar to riboflavin.

Bioinformatics software and databases are essential components of modern life science research, yet their mentions in the scientific literature are often inconsistent and difficult to systematically identify at scale. The lack of a comprehensive and up-to-date catalog of bioinformatics resources hinders efforts toward automated biomedical knowledge extraction and streamlined data analysis. Here we present SNAIL, a hybrid named entity recognition framework designed to automatically identify bioinformatics software and database (SW/DB) names from biomedical texts. SNAIL integrates complementary lexical and semantic modeling strategies. The lexical component captures orthographic patterns and contextual cues characteristic of SW/DB names, while the semantic component leverages contextual embeddings generated by transformer-based language models such as SciBERT, combined with an explicit token-masking strategy to enhance entity-focused representations. A large training corpus was constructed automatically through a hybrid pipeline that integrates citation-hinted extraction with large language model-assisted distillation. Evaluation on two independent benchmark datasets and real-world research articles demonstrates that SNAIL substantially outperforms existing approaches, including domain-specific methods such as bioNerDS2 and general-purpose large language models such as ChatGPT, Gemini, Grok and Claude. Applying SNAIL to large-scale literature analysis further reveals distinct journal-level preferences across bioinformatics subfields. These results demonstrate that SNAIL provides an accurate and scalable solution for identifying bioinformatics resources in scientific texts and enables systematic meta-analysis of tool usage and research trends.

Background: Reproducible taxonomic collapsing and geological-timescale annotation of time-calibrated phylogenetic trees in R often require coordination among several packages and repeated code for label parsing, clade validation, plotting, and export. Workflow-managed analyses additionally benefit from non-interactive configuration, predictable diagnostics, and machine-readable exit status. Results: We present Rclade, an R package that consolidates the multi-package coordination required for taxonomic collapsing into a streamlined, single-function interface. Rclade provides (1) custom ggproto objects (GeomPolygonStraight/GeomSegmentStraight) that bypass coord_munch() interpolation to achieve straight-edge rendering of collapsed triangles in circular layouts; (2) automatic detection and parsing of four taxonomic-label formats (GTDB, Silva, NCBI, embedded) plus user-supplied custom regex, with explicit input-validation contracts and parsing-accuracy evaluation on real and derived test sets; and (3) workflow embeddability through YAML configuration, library-mode APIs, and standard Unix exit codes. Benchmarks on synthetic and real datasets (200-10,000 synthetic tips and real reference trees up to 10,122 tips; 5 replicates at every scale under a unified fully rendered measurement protocol) show that the full-pipeline overhead is modest for interactive use (median {approx}0.87 s in-session rendering and {approx}8.4 s process-level wall-clock at 10,000 tips). Conclusions: Rclade is a convenience layer over the ggtree/deeptime ecosystem that reduces boilerplate while adding targeted technical improvements for circular-layout rendering and format heterogeneity management.